arXiv: 1504.0761 lv2 [hep-th] 15 May 2015 


MIT-CTP /4668 


May 18, 2015 

Off-shell hydrodynamics from holography 

Michael Crossley, Paolo Glorioso, Hong Liu and Yifan Wang 
Center for Theoretical Physics, 

Massachusetts Institute of Technology, 

Cambridge, MA 02139 

Abstract 

We outline a program for obtaining an action principle for dissipative fluid dynamics by considering 
the holographic Wilsonian renormalization group applied to systems with a gravity dual. As a 
first step, in this paper we restrict to systems with a non-dissipative horizon. By integrating out 
gapped degrees of freedom in the bulk gravitational system between an asymptotic boundary and 
a horizon, we are led to a formulation of hydrodynamics where the dynamical variables are not 
standard velocity and temperature fields, but the relative embedding of the boundary and horizon 
hyper surfaces. At zeroth order, this action reduces to that proposed by Dubovsky et al. as an 
off-shell formulation of ideal fluid dynamics. 
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I. INTRODUCTION 


At distance and t im e scales much larger than the inverse temperature and any other 
microscopic dynamical scales, a quantum many-body system in local thermal equilibrium 
should be described by hydrodynamics. Except for ideal fluids, the current formulation of 
hydrodynamics has been on the level of equations of motion. There are, however, many 
physical situations where hydrodynamical fluctuations play an important role. An action 
principle is greatly desired. There are two main difficulties. One is to properly treat dissi¬ 
pation, and the other is to find the right set of dynamical degrees of freedom to formulate 
an action principle, as standard variables such as the velocity field appear not suitable. 

In principle it should be possible to derive hydrodynamics as a low energy effective held 
theory from a quantum held theory at a finite temperature via Wilsonian renormalization 
group (RG) by integrating out all gapped modes, but in practice it has not been possible to 
do so. Such a formulation should lead to an action principle. 

For holographic systems, the holographic duality [1-3] provides a striking geometric de¬ 
scription of the renormalization group how in terms of the radial how in the bulk geometry. 
In particular, the Wilsonian renormalization group how of a boundary system can be de¬ 
scribed on the gravity side by integrating out part of the bulk spacetime along the radial 
direction [4, 5] . The proposal expresses the Wilsonian effective action in terms of a gravita¬ 
tional action defined at the boundary of the remaining spacetime region. 

The purpose of the current paper is to take a hrst step toward deriving an action for 
hydrodynamics using holographic Wilsonian RG. 1 The basic idea is as follows: consider the 
gravity path integral 

Z[g 1} g 2 }= f DGe iS ^ (1.1) 

JM C 

over metrics G on a complex bulk manifold A4 C consisting of two copies of asymptotic AdS 
spacetimes patched together at a dynamical horizon hypersurface, as shown in Fig. 1 . A4 C has 
two asymptotic boundaries i hM 1,2 with boundary metrics gi,g 2 respectively. The horizon 
1 The connections between holography and hydrodynamics has by now quite a history, starting with [6-9] 
and culminated in the fluid/gravity correspondence [10-12], 
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FIG. 1. (a) Complex bulk manifold A4 C consisting of two copies of asymptotic AdS spacetimes 
patched together at a horizon hypersurface. Also labeled are stretched horizons £i ,£2 discussed 
around (1.2). (b) A boundary theory Schwinger-Keldysh contour used to describe non-equilibrium 
physics. The two AdS regions map to the two horizontal legs of the Schwinger-Keldysh contour, 
while the analytic continuation around the horizon corresponds to the vertical leg which defines 
the initial thermal density matrix. 


is dynamical as its metric is integrated over. The two copies of AdS can be considered 
as corresponding to the two long horizontal legs of a Schwinger-Keldysh contour, with the 
continuation around the horizon corresponding to the vertical leg [13]. In (1.1) one integrates 
out all gapped degrees of freedom, and the resulting effective action for whatever gapless 
degrees of freedom remain is the desired action for hydrodynamics. For this purpose, it is 
convenient to introduce stretched horizons £„,« = 1,2 on each slice of the bulk manifold, 
which separate the bulk manifold into three different regions (see Fig. 1), i.e. 

/•Si rYi 2 rdMc 

I + / + / • (l- 2 ) 

dM\ Jt, i Jt , 2 
The bulk path integral can be written as 



Z{gi,g2\— / Dh\ / Dh 2 4'ir[^i, h 2 \ Tuv[/ii, (h] v biuv[h2, # 2 ] 
J Si J S 2 


(1.3) 


where h\ and h 2 are induced metric on the stretched horizons. Various factors in the inte¬ 
grand of (1.3) arise from the path integrals in three regions, e.g. 

rhi 


*uv[hi,9i]= / DGe 13 ^ 


(1.4) 


'9 1 
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integrates over all metrics G between dM.\ and Si with Dirichlet boundary conditions g\ 
and hi, and similarly with the others. The complex conjugate on 'hu V [h 2 ,^ 2 ] is due to that 
the bulk manifold in the region between S 2 and dAi 2 has the opposite orientation from that 
between dAi l c and Si. 

Connections between hydrodynamics and Schwinger-Keldysh contour have been made 
recently in various contexts in [14-21], 



FIG. 2. The gapless degrees of freedom in the path integrals (1.4) are the relative embedding 
coordinates Xf(cC) between the horizon and the boundary. Xf can be understood geometrically 
as follows: start with cC at Si, shooting a congruence of geodesics orthogonal to Si toward the 
boundary, the intersections of these geodesics with the boundary define Xf. 

In this paper we describe integrations over gapped degrees of freedom in the path inte¬ 
gral (1.4). As anticipated earlier by Nickel and Son [22], in (1.4) the only gapless degrees 
of freedom are the relative embedding coordinates A"“(cr M ) of the boundary Ai l c and the 
stretched horizon hypersurface Si, see Fig. 2. Integrating out all other degrees of freedom 
we obtain an effective action iu V [Xf, gi, hi] for embeddings X“, i.e. (1.4) becomes 

Vuv[hi,gi) = I DX?e iIvv[x ?’ Sl ’~ hl] . (1.5) 

We develop techniques to compute /uv \X“,gi, hi\ in expansion of boundary derivatives at 
full nonlinear level in a saddle point approximation. Plugging (1.5) into (1.3) and evaluating 
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hi, Ii 2 integrals one then obtains the full hydrodynaniical action in terms of and X $, i.e. 

Z[gi,g 2 \ = J DX“DX2e iIh y d '° [Xl ’ Xr ’ Sl ’ S2] (1.6) 

with 

e i/hyd»[x 1 ,x 2 ;g 1 , S 2] = f J D^ 2 ^ IR [^ 1) ^ 2 ]e i ( J uvAi Q >9i:H--fuv[x 2 “,ff 2 ,/ l2 ]) _ (i.7) 

J Si J S 2 

The evaluation of Tir [hi, J 12 ] requires developing new techniques for analytic continuations 
through the horizon. We will leave its discussion and the full evaluation of (1.7) elsewhere. 
Hydrodynaniical actions based on doubled X a degrees of freedom discussed here have also 
been discussed recently in [16-21], 

We also show that at zeroth order in the derivative expansion, if one (i) takes h\ to the 
horizon, i.e. making Ej a null hypersurface, and (ii) requires hi to be non-dissipative, i.e. 
the local area element is constant along the null geodesics which generate the horizon, h\ 
completely decouples from /uv[AT“, gi, hi], and Juv reduces to the conformal version of the 
ideal fluid action proposed by Dubovsky et al [23, 24], i.e. 

/uvK, <7l) hi] = /ideal[£l, 9l] (1-8) 

where 

4ieai [£,<?] = -(d- 1) J d d ay/—^ (det cW 1 ) , (« _1 ) 1J = d^d^, (1.9) 

and £‘((7^) (with i = 1, 2, ■ • - , d — 1) are embeddings X a (cr M ) for a null hypersurface for 
which the time direction decouples. In particular, the volume-preserving diffeomorphisms 
which played a key role in the formulation of [23] arise here as residual freedom of horizon 
diffeomorphism. The entropy current also arises naturally as the Hodge dual of the pull-back 
of the horizon area form to the boundary. 

It is tempting to ask whether conditions (i) and (ii) in the previous paragraph will also 
lead to a non-dissipative fluid action at higher orders. We find, however, that the 2nd order 
action is divergent unless one is restricted to shear-free flows. While it makes sense to make 
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such restrictions in an equation of motion, imposing it at the level of path integrals for X a 
appears problematic. We thus conclude that one must include dissipation in order to have 
a consistent formulation. 

We also note that the fact we find (1.9) when pushing Si to the horizon does not neces¬ 
sarily imply that at zeroth order the full effective action (1.7) will be given by 

7}iyd ro /ideal [7l • fj ]\ 7 deal [^2* i/d] (1.10) 

as the integrations over hi, I 12 will generate new structures. At this stage the precise relation 
between (1.9) and the zeroth order of /hydro is not yet fully clear to us. 

The plan of the paper is as follows. We will explain our holographic setup and the 
gravitational boundary value problem in details in Sec. II. In Sec. Ill, we perform the path 
integral (1.4) using saddle point approximation to obtain Juv defined in (1.5) at zeroth order 
in the boundary derivative expansion and relate it to the ideal fluid action (1.9). In Sec. IV 
we briefly comment on the generalization to higher orders in the derivative expansion. We 
conclude with a discussion of open questions and future directions in Sec. V. 

Note added: We understand Jan de Boer, Natalia Pinzani Fokeeva and Michal Heller have 
obtained similar results [25]. 


II. SETUP 

In this section, we describe in detail our setup for computing (1.4) to obtain / uv [AT a , g, h] 
defined in (1.5). 


A. Isolating hydrodynamical degrees of freedom 

I 11 this subsection, we describe a series of formal manipulations of path integrals for 
gravity which allow us to isolate X a as the “would-be” hydrodynamical degrees of freedom. 
There are standard difficulties in defining rigorously path integrals for gravity, which will 
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not concern us as we will be only interested in the path integrals at a semi-classical level, 
i.e. in terms of saddle points and fluctuations around them. 

Consider a path integral of the form 

fh 

V[h,g]= / DGe lS[G] (2.1) 

where the integration is over all spacetime metrics 

ds 2 = G ! Mv( cr )dcr A/ dcr A ' = N 2 dz 2 + + N tl dz)(d(j u + N u dz) (2.2) 

between two hypersurfaces Suv and Sir at some constant-^; slices and whose respective 
intrinsic geometries are specified by g /w and h^, i.e. 

X/iv | Suy )j — ^7-“ / ( cr ) • (2-3) 

In (2.1), one should integrate over all values of N and iV^ without any boundary conditions 
for them on Suv and Sir. For this paper we will only be concerned with evaluating (2.1) to 
leading order in the saddle point approximation, thus will not be careful about the precise 
definition of integration measure, ghosts, and Jacobian factors for changes of variables. We 
will comment on these issues in Sec. V. 

Later we will take Suv to the boundary of an asymptotic AdS spacetime and Sir to an 
event horizon. We will for now keep them arbitrary for notational convenience. We will 
also for now keep the gravitational action S'fG] general assuming only that it is diffeomor- 
phism invariant and that the boundary conditions (2.3) give rise to a well defined variational 
problem. The variation of the action then has the form 

6S= l j dd+l(T ^-^ EMN8G MN (2.4) 

without any boundary term. The equations of motion are thus 

E mn = 0 (2.5) 

while diffeomorphism invariance of the action S[G] also leads to the Bianchi identities 

MN _ q 
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VmE 


( 2 . 6 ) 



Given the diffeo morphism invariance of the bulk action S[G], ty[h,g] is invariant under 
independent coordinate transformations DiffjR x Diffuv of the two hypersurfaces, i.e. 

^[h,g] = ^/[Aj R /iA ir , Auy^A uv ] (2.7) 

where Air and Auv denote independent coordinate transformation matrices on Eir and Euv- 
Now consider transforming the metric (2.2) to the Gaussian normal coordinates (GNC) 
£ a (ct) = (u,x a ) in terms of which 

ds 2 = du 2 + 'y a b{u, x)dx a dx b = GABd^ A d^ B . (2.8) 

Here we choose Gaussian normal coordinates for later convenience. The subsequent dis¬ 
cussion applies with little changes to any set of “gauge fixed” coordinates. The metric 
components Gmn can be expressed in terms of ( 7 ^, £ A ) as 

Gmn = Gab9m^ A 9n^ B = OmuOnu + ^ ab dMX a dN'X b ■ (2.9) 

In choosing the Gaussian normal coordinates we have the freedom of fixing the values of u 
and x a at one end. For our later purpose it is convenient to choose a hybrid fixing 

uL = Mo = const, x“L = (2.10) 

I ^uv 1 ^IR H 1 V 7 

The values of u at Sir and x a at Euv are then determined dynamically, which we will 
parameterize as 

« 4 m = f K)> He„ v = . (2.11) 

In terms of the foliation of (2.8), Euv and Eir can thus be written as 

S uv : u = Uqj Sir : m = r(x“) = r(X _ 1 (a: a )) ( 2 . 12 ) 

and the boundary conditions (2.3) now become 

b ab \ u =r(x a ) = ^ a6 ’ Ta 6 1 U=UQ = dab (2-13) 

with 

r)(j v f)X a 

hab = Kb - d a rd b T , g ab (X) = — — = {r lt gj ~ l ) ab , J% = — . (2.14) 
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Note that = (u(a AI ),x a (cr AI )) are dynamical variables and in going from (2.2) to (2.8), 
we have essentially traded Gmn = (N, N^Xnv) f° r (w, £ a , 7 a&)- The path integral can now 
be written as 


V[h,g] = Dx a Du D lab e 


iS[u,x a ,"/ ab \ 


(2.15) 


where 7 is required to satisfy the boundary conditions (2.13). The coordinate invariance of 
the action implies that the action is independent of the bulk fluctuations of u and x a . Thus 
the path integrals over x a and u reduce to those over the boundary fluctuations ( 2 . 11 ) 

r r ph 


Mh,9}= / DX° / Dr / D lab e 


iS[r, 7] 


(2.16) 


In the above integrals X a always appears with g through the induced metric g defined 
in (2.14). In addition to appearing in the IR boundary condition h ab for ~f ab integrals, r also 
appears in the action S explicitly as the IR integration limit and boundary terms (which we 
will specify below). 

X a and r can be considered as the “Wilson line” degrees of freedom associated with 1V M 
and N. Physically X a (a a ) describes the relative embedding between the coordinates x a on 
£ ir and the coordinates cr M on Euv, while t(x) describes the proper distance between S uv 
and S IR . 

The path integrals (2.16) will be evaluated in stages: we first integrate over all possible 
7 ab with a fixed relative embedding X a and proper distance r to find 


Jh[r,h,g] _ 


= / Djah e 


iS[r,7] 


(2.17) 


and then integrate over r (i.e. all possible proper distances) 


e iluv[h,g] 


J Dr e iIl[rX9] . 


(2.18) 


In the path integrals (2.17), with a finite r, u direction is essentially compact and 7 a b can be 
consistently integrated out to yield a local action Ii [r, h, g ], which can be expanded in the 
number of boundary derivatives of r, h and g , assuming they are slowly varying functions 
of boundary coordinates. In the boundary theory language 7 a b should thus correspond to 
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modes with a mass gap. r depends only on boundary coordinates, and does not contain 
derivatives at leading order, and thus can also be consistently integrated out. By definition 
X a always come with boundary derivatives as in (2.14), i.e. they correspond to gapless 
modes, and thus should be kept in the low energy theory. Integrating them out will lead to 
nonlocal expressions. 

Let us briefly consider the symmetries of Ii[r,g,h]. It is invariant under ^-independent 
diffeomorphisms of x a —> x' a (x) under which g, h transform simultaneously as tensors. These 
are large diffeomorphisms as they change the asymptotic behavior of AdS. I\ is also invariant 
under diffeomorphisms of as it contains X a and g only through g, which is invariant due 
to the canceling of transformations in g and X a , 


9rW) 9»M) = -^Ap(^V)) q-;, X a (a) = X a (a\a)) . (2.19) 


This implies that 



boundary stress tensor, we then see that the X a equations of motion are equivalent to the 
conservation of the boundary stress tensor. Parallel statements can be made about Juv 
which comes from integrating out r. 

Finally to conclude this subsection, let us be more explicit about the UV boundary 


condition. In an asymptotic AdS spacetime with Euv at a cutoff surface near the boundary, 
7 a b in ( 2 . 8 ) should have the behavior 



( 2 . 21 ) 


with L the AdS radius and r y < ^(x a ) finite. Thus we should replace the boundary condi¬ 
tion (2.13) at u = mq by 



( 2 . 22 ) 


where 


9ab — ( J gJ ) ab 


(2.23) 
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and g MI ,(cr) is the background metric for the boundary theory. 


B. Saddle point evaluation 


Now consider evaluating the path integrals (2.16)~(2.18) using the saddle point approxi¬ 
mation. To elucidate the structure of equations of motion for r and X a , we consider (2.4) 
now with Gmn considered as a function of ^ ab and £ A via (2.9). Under variations of 7^, we 
have 

5G M n(x) = S~/ ab d M x a d N x b (2.24) 

which then (2.4) implies the equations of motion 

E MN d M x a d N x b = 0, => E ab = 0 (2.25) 


where E ab = 0 is the aft-component of the equations of motion in coordinates (2.8). Below 
we will refer to these equations as “dynamical equations.” 

Under variations of £ A , we have 

5Gmn = q^c ^ 4 B S^ C + 2 GAcdMi A dN^ C ■ (2.26) 


The bulk part of (2.4) then leads to the Bianchi identities in coordinates (2.8), which is as 
it should be since £ a (ct) is a coordinate transformation. But now there are boundary terms 
remaining 


dz 


5S — I d d xV^HE A B ^8£ A 


<9£ 


dz 


Sir 


d d xyf^HE A B —5t A 


(2.27) 


Suv 


which upon using (2.10)-(2.12) implies that 

dz 


y/—HE u B 


<9£ 


B 


= 0 =* \E U U - E u a 


dr_ 

dx a 


= 0 


(2.28) 


and 


V^HE a B 


dz 


dt 


B 


= 0 


E a U L =0 

I^UV 


(2.29) 
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with (2.28) corresponding to the equation of motion from varying r while (2.29) corresponds 
to those from varying X a . In deriving the second equations in both (2.28) and (2.29) we have 
assumed that \J—H and at Sir and Suy are nonzero. It can be readily checked that (2.25) 
and (2.28)-(2.29) are equivalent to (2.5), and that the Bianchi identity ensures that (2.28) 
and (2.29) are satisfied everywhere once they are imposed at Sir and Suv, respectively. 
Following standard convention, below we will refer to (2.28) as the Hamiltonian constraint 
and (2.29) as the momentum constraints. 

Now recall from the general results on the holographic stress tensor [26] that the momen¬ 
tum constraints (2.29) in fact correspond to the conservation of the boundary stress tensor 


V a T ab = 0 


(2.30) 


where V a is the covariant derivative associated with g ab and T ab is the stress tensor for the 
boundary theory with background metric g ab in the state described by (2.8). Since g )W and 
g a b a re related by a coordinate transformation equation (2.30) is equivalent to 


= 0 


(2.31) 


where is the covariant derivative associated with g^ and is the stress tensor for the 
boundary theory with background metric g^. This gives an alternative way to see that X a 
equations of motion are equivalent to conservation of the boundary stress tensor. 

At the level of saddle point approximation, A as defined in (2.17) is obtained by solv¬ 
ing (2.25) for 7 ab , and Juv in (2-18) by solving (2.28) for r(x a ). In other words, Juv is 
computed by evaluating the gravity action with dynamical equations and the Hamiltonian 
constraint imposed, but not the momentum constraints. 
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C. Einstein gravity 


We now specialize to Einstein gravity, in which the gravity action A[G] in (2.1) can be 
written in the Gaussian normal coordinates (2.8) as 

S[T, 7 ] = f d d+1 i y/=j(R - 2A) - [ 2K + [ 2 K + S ct 

J J Sttv J Etr 


d d+1 ti a/— 7 [ {d) m - (K ab K ab - A' 2 ) - 2A] + S m + A ct 


(2.32) 


where 


K ab = 1 -d nl « b = \i al> , A-idyVn. K = K\= 1 -r b lU (2.33) 

are extrinsic curvatures for a constant-w hypersurface, and S ct is the standard AdS countert¬ 
erm action at the Erv [26] 

5c, = / <fV=7 (+ •••)■ (2.34) 

J u=UQ^—oo \ Jj CL Z J 

Air is a boundary action at S IR which arises from the fact that E IR , given by u = r(x a ), 
is not compatible with the foliation of constant-?/ hypersurfaces, and can be written as (see 
Appendix A for a derivation) 

1 
n 


Sir — 2 


with 


d d xV~h - (( h ab (dr ) 2 - d a Td b r) K ab - D 2 t ) 
n = \/l — d a rd a T 


(2.35) 


(2.36) 


where the indices are raised and lowered by the intrinsic metric h ab on Sir and D is the 
covariant derivative associated with h ab . 

For convenience below we will use /C to denote the matrix K a b and thus K = Tr/C. 
Various components of the Einstein equations in Gaussian normal coordinates (2.8) can 
then be written as 

-El = K! - Tr/C' + /CTr/C - Ric (d) [ 7 ] - ^ (Tr/C 2 + Tr 2 /C) + ^ R {d) [ 7 ] - A = 0 (2.37) 

_ E uu = I Tr /c 2 _ ^Tr 2 /C + ^R {d) [ 7 ] - A = 0 (2.38) 

~E u a = VJ< - V b K b a = 0 (2.39) 
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with V a the covariant derivative associated with 7 ^. As discussed in Sec. IIB, in order to 
not impose conservation of the stress tensor, i.e. leave hydrodynamical modes off-shell, at 
the saddle point level we should not impose the momentum constraint (2.39). We only need 
to solve the dynamical equations (2.37) for 7 ab and a combination of (2.38)-(2.39) at Sir for 
r (see (2.28)). 

From now on we will set the AdS radius L — 1. 

III. ACTION FOR AN IDEAL FLUID 

In this section we first evaluate explicitly I\jy[X a ,g, h} defined in (1.5) at zeroth order in 
the derivative expansion, assuming that X a , g, h are slowly varying functions. We then show 
that pushing h to a horizon hypersurface and requiring it to be non-dissipative, we obtain 
the ideal fluid action of [23, 24], 

A. Solving the dynamical equations 

We will perform the 7 integrals (2.17) using the saddle point approximation, i.e. it 
boils down to solving the dynamical Einstein equations (2.25) at zeroth order in boundary 
derivatives. At this order we can neglect boundary derivatives of t(x), J a M , and 7 The 
boundary conditions for 7 ab become 

7 (u = r)=h = h, 7 (tt = u 0 -A - 00 ) = e~ 2u °g, g ab = (j _ 1 ^J _1 ) ab . (3.1) 

For notational simplicity here and below we will often use 7 and g to denote the whole matrix 
7 ab and g ab . Hopefully the context is sufficiently clear that they will not be confused with 
their respective determinants. 

Explicit expressions for the Einstein equations in Gaussian normal coordinates (2.8) are 
given in Sec IIC. At zeroth order in boundary derivatives, the dynamical part of the Einstein 
equations (2.25) (more explicitly (2.37)) becomes 

K! - Tr/C' + /CTr/C - ^ (Tr/C 2 + Tr 2 /C) - A = 0 (3.2) 
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which can be rewritten as 


-—- (K' + - = -A - - TrK 2 , 

d V 2 ) 2 


K' + KK — 0 

where IK is the traceless part of /C 


From (3.4) 


/C = K + — 1, TrIK = 0 
d 


-(TrK 2 )' = -A'TrK 2 


Taking derivative on both sides of (3.3) leads to 

d- 1 


d 


(K" + KK’) = A'TrK 2 


(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


Eliminating TrK 2 between (3.6)-(3.7) and using (3.3) we then find an equation for K 

K" + 3K'K + K(K 2 -d 2 ) = 0 (3.8) 


which is solved by 


K = d- 


0l\Cl2£ 


2 du ^ 


(3.9) 


(1 + Oiie du )(l + a 2 e du ) 

where are some constants. Inserting (3.9) into (3.4) we find 

K _ 5 0 (aq - a 2 )e du 

(1 + aie du )(l + a 2 e du ) 

with b 0 a constant traceless matrix. Plugging (3.10) into (3.7) we find that b 0 has to satisfy 


(3.10) 


Tr6g = rf(rf — 1) . 


Combining (3.9) and (3.10) we then obtain 

1 


K. = 


(&o(oa — cx 2 )e du + {ot\OL 2 e 2du — 1)) 


(3.11) 


(3.12) 


(1 + aie du )(l + a 2 e du ) 

Now integrating (2.33) and imposing the boundary condition at UV (i.e u = u 0 —* —oo) 
we find that 

7 = ge~ 2u ((1 + aie d “)(l + a 2 e du )) 1 e 2k(u)b ° (3.13) 
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with 

i , \ 1 , (1 + a i eC 

k(u) = - log --, 

d \ 1 + <y 2 e du 

where here and below we will always take aq > c^- 
Introducing 

H = g _ 1 /z, oi\ c = aqe dr , a 2c = ol 2 e dr 
from the IR boundary condition 7 ( 7 / = r) —h we then hnd 


(3.14) 


(3.15) 


Vdet// - e dr (l + Qi c )(l + a 2c ), 6 0 - 777 , log Id, Tr(logld) 2 - — , ^ log 2 1 + Ql 


2k(r) 


d 


1 + ol 2c 
(3.16) 


where H denotes the unit determinant part of H and the last equation of (3.16) follows 
from (3.11). Requiring the metric 7 a f, to be regular and non-degenerate between u = — oo 
and r, we need 

1 + ol\ c > 0, l + a 2 c >0. (3.17) 

Given H and r, we can use the first and last equations of (3.16) to determine a\ )2 and 
then the second equation to find b 0 . Note at this stage r is not constrained by H and thus 
can be chosen independent of h. More explicitly, 


«i c = (det H) 4 e 2 e 2 —1, a 2c = (det H ) 4 e 2 e 2 — 1 


with 


Zr - 


d 


4(d — 1) 


Tr(log H'y 


(3.18) 


(3.19) 


B. Effective action for r and X a 


At zeroth order in boundary derivatives, the Einstein action (2.32) becomes 
S[t, 7 ] = I d d+1 £ a /— 7 [- (Tr K 2 - K 2 ) - 2A] + 5 cl 
and substituting (3.12) into (3.20) we have 

Ii [r, h , 51] = 2 (d - 1) f d d x\f zr g [-e~ dr + e _cZA + aia 2 e dr ] + 


(3.20) 


(3.21) 
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with the counterterm action given by 


S ct = -2 (d - 1) / d d x^ = -2 (d - 1) 

We then find that 

h[r, h,g\ = -2(d - 1) J d d xV zr g Ci(H,r) 

with 


J d d xV~g (e dA —a.\ — a 2 ) 


(3.22) 


(3.23) 


Hi ( H , t) — e -cvi -1- c^2 — oticx 2 6 

= —VdetH + 4e - 5 dT (det H)* cosh ^ — 2e“ dT (3.24) 

where in the second line we have expressed the integration constants a in terms of bound¬ 
ary conditions via (3.18). We notice that at zeroth order, I\ depends on h and g only through 
the combination H = g~ 1 h. This follows from that I\ must be invariant under the diffeo- 
morphisms of x a for which h and g transform simultaneously, as noted in the paragraph 
before (2.19). At zeroth order Tr H n for n = 1, 2,..., d are the only independent invariants, 
r can now be integrated out by extremizing A which gives 

e~ dT ° = VdetH cosh 2 y (3.25) 


and thus 

£uv = £i(H,t 0 ) = VdetH cosh^ c . (3.26) 

Collecting everything together we thus find that 2 

1 

(3.27) 


Juv[A" a ; h, g] = —2 (d — 1) J d d xV~g V det H cosh jj Tr(log H)‘ 


One can readily check that the same result is obtained by solving instead the Hamiltonian 
constraint (2.28) at zeroth order. Also note that with r = r 0 given by (3.25), equation (3.18) 
becomes 


oi\ c — — oi 2c — tanh 


(A 

2 


(3.28) 


2 The overall minus sign has to do with our choice of orientation of bulk manifold. 
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C. Horizon limit 


Equation (3.28) implies that after integrating out r, a 2 = — cq is negative. We will now 
simply rename «i as a. It is convenient to introduce 


a = e dUh 


(3.29) 


then from (3.28) 


1 . . z c 

u h = r - - log tanh — > r . 
a 2 


(3.30) 


Now if we extrapolate the solution (3.13) beyond u = r all the way to Uh, then ^/— det 7 
develops a simple zero at u = Uh, which we will refer to as a “horizon.” Since the “horizon” 
lies outside the region where our Dirichlet problem is defined, 7 a b does not have to be regular 
there, so this does not have to be the horizon in the standard sense. Now let us consider a 
sequence of h whose time-like eigenvalue approaches zero. Equivalently, an eigenvalue of H 
(which we will denote as h 0 ) goes to zero. At h 0 = 0 , h describes a null hypersurface and 
Sir becomes a horizon for the metric between Sir and Euv- We thus define the h 0 —y 0 
limit as the hydrodynamic limit. 

In this limit, we have 


det H —> hnSdetH , z c —> — log h 0 H——-- 

0 ’ c 2 & 0 2(d — 1 ) 


log Sdet H 


(3.31) 


where Sdet H denotes the non-vanishing subdeterminant of H and can be written as 


Sdet H = Vd _ x (Tr H, Tr H 2 , ■ ■ ■) 


(3.32) 


where p d ~i is the standard polynomial which expresses the determinant of a non-singular 
(d — 1) x {d — 1) matrix in terms of its trace monomials. From (3.25) and (3.29) we thus 


End 


e dT ° —> -(Sdetih) 2 ( d - 1 >, Uh — r —>■ 0 

and the action (3.27) becomes 

(3.33) 

/uv — — (d — 1 ) j 

f d d x\f^g (Sdet H) 2 ( d - 1 ) 

(3.34) 

T-1 

1 

T 

f d d cr (Sdet/7) 2 ( d -!) . 

(3.35) 
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D. Entropy current 


Here we discuss the geometric meaning of Sdet H and (3.35). Denote the null eigenvector 
of h by £ a , which give rises to a congruence of null geodesics which generate the null hyper¬ 
surface. We can then choose a set of coordinates ( v , p) on Sir with v the parameter along 
the null geodesics generated by i a and £* remaining constant along geodesics. In this basis, 
we then write the metric on Sir as 

Qp Qp 

ds| IR = habdx a dx b = cry (v,£)dpdp, h ab = i = 1,2, ■ ■ - ra, n = d- 1 . (3.36) 

It then follows that 

Sdet H = det(a lj crjk) = det a det oT x (3.37) 

where ct -1 is dehned as 

Qp ftp 

(a-'r = ^ ^ J\da- = ^ . (3.38) 

We thus find that VSdet Ed can be written as horizon area density y/o normalized by the 
“area density” of the pull back of boundary metric g to Sir. 

The physical meaning of r/Sdet H can be better elucidated if instead we pull back all 
quantities to the boundary. We now show that it can be interpreted as a definition of (non¬ 
equilibrium) entropy density of the boundary system. 3 For this purpose, consider the area 
form on the Sir which can be written as 

a = y/a d£ ] A d£ 2 A • • • A dp 1 . (3.39) 

Note that the horizon area ydx has no physical meaning itself in the boundary theory as 
its definition depends on a choice of local basis. It does become a physically meaningful 
quantity when we pull it back to the boundary via the relative embedding map J a M introduced 
in (2.14). More explicitly, 

a = yfdE 1 A E 2 A • • • A E n , E 1 = E\da^ . (3.40) 

3 The construction below is similar to the construction of entropy current in fluid/gravity [12, 27]. 
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From (3.40) we can define a current which is the Hodge dual of a on the boundary 


jV — 6 W" v n ( 


^1 —Vn 


- . rph ... pi 


TV. 


(3.41) 


where e ,w ''" Vri is the full antisymmetric tensor for g and e^...^ is that for cry. Similarly, it is 
natural to pull back the null vector £ a to the boundary, giving 


= (J-y a f\ g^uV = -1 


(3.42) 


and we have chosen a convenient normalization for u R By construction, is parallel to 
and we can then write 

f = su ", s 2 = ■ (3.43) 

From (3.41), we hnd that 

s = \J det(a 1 jajk) = VSdetH . (3.44) 

We will interpret as the velocity field of the boundary theory, j^ (divided by 4 Gn) as the 
entropy current, and s (divided by dGjy) as the local entropy density. All these quantities 
are independent of choice of local coordinates on Sir. We also stress that their definitions 
do not depend on the derivative expansion and thus should apply to all orders. 

With this understanding the action (3.35) can be written as 

Juv = ~ J d d a\f-~H e(s) . (3.45) 

where 

e(s) = (d — l)^ 3111 (3.46) 

has precisely the scaling of the local energy density as a function of entropy density for a 
conformal theory. From the perspective of evaluating the bulk action it can also be under¬ 
stood as follows: the bulk integration in (2.32) can be interpreted as giving the free energy 
while the Gibbons-Hawking term at the IR hypersurface becomes equal to entropy times 
temperature in the horizon limit. Their sum then gives the energy of the system. 
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E. Hydrodynamical action and volume-preserving diffeomorphism 


We now impose that the system is non-dissipative, which amounts to requiring that the 
entropy current (3.41) is conserved 

= 0. (3.47) 

where V is the covariant derivative for g. The above equation can also be written equivalently 
in various different ways in terms of horizon quantities. In terms of the horizon area form 
(3.39), 

da = 0 (3.48) 

or area density 

<9„det<7 = 0 —>■ det cr = det cr(£) . (3.49) 

i.e. the horizon area is independent of the horizon “time” v. Note that the form of the 
metric (3.36) is preserved with a u-independent coordinate transformation 

r e"(0 (3.50) 

which can be used to set 

deta = l —» SdetH = det a -1 . (3.51) 

We thus see that with non-dissipative boundary condition at zeroth order the horizon metric 
completely decouples in the hydrodynamical action, and we find 

Juv = -(d- 1) J d d ay/—§ (detcT 1 )^ , (oT 1 )^ = (3-52) 

which is precisely that of [23, 24] applied to a conformal theory. 

After fixing (3.51), there are still residual volume-preserving diffeomorphisms in which 
played an important role in the formulation of [23, 24]. Here they arise out of a subgroup 
of horizon diffeomorphisms which leave “gauge fixing condition” (3.51) and the coordinate 
frame (3.36) invariant. If the non-dissipative horizon condition (3.49) can be consistently 
imposed to higher orders in derivative expansion, we should expect the resulting higher 
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order action to respect the volume-preserving diffeomorphisms. As we will see in Sec. IV, 
however, at second order we encounter divergences, which implies that (3.49) can no longer 
be consistently imposed for unconstrained integrations of X a . 


F. More on the off-shell gravity solution 


Here we elaborate a bit further on the off-shell gravity solution (3.13). First let us collect 
various earlier expressions. After integrating out r, with (3.28) and (3.29) the off-shell 
metric (3.13) can be written as 

7 = ge~ 2u (1 - a 2 e 2du Y e^ log ^ (3.53) 


where we have introduced 


z(u) = log 


1 + ae 


du 


1 — ae' 


du 


—logtanh -(«/, — u) = tanh ae' 


du 


(3.54) 


Also recall that 


e dT = y/( det H cosh 2 -A, a — e dUh = - \J det H sinh z Cl A = d(uh — r) = — logtanh -A . 

(3.55) 

With a general regular h, the solution when extrapolated to the “horizon” Uh is singular. 
This is perfectly okay as physically the behavior of the metric outside the region is of no 
concern to us. We will now show that if we do require the extrapolated metric to be also 
regular at the “horizon” u — Uh, i.e. u = Uh becomes a genuine horizon, then we recover the 
standard black brane solution. Of course this also implies that h has to take a very specific 
form. 

We now impose a “regularity” condition: has only one eigenvalue approaching zero as 

the horizon is approached with the other eigenvalues finite. Near the horizon, 5 = d(uh~u ) —y 
0 with 

£ 

z(u ) — y — log - —> +oo . (3.56) 
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Denoting the eigenvalues of log H and g *7 as b M and 7 ^ respectively, from (3.53) we then 
have 

_ bj± 

(0 “ . (3.57) 

If we denote the time-like eigenvalue of g~ l 7 by 70 and the rest by 7 $, the regularity condition 
amounts to that 70 goes to zero while all 7 * finite. Requiring 7 \ to be finite implies that 

k = = -^logtanh^ = b, b 0 = = -{d - l)b (3.58) 

where the second equation follows from that log H is traceless. We thus find that the system 
has to be isotropic! 

Denoting the time-like eigenvector vector of H as £ a , from (3.58) we can write H as 

H\ = (d - 1 )U a i b + b(5 a b + t£ b ) (3.59) 

where we have defined 

4 = g ab £ b , £ a L = -1 • (3.60) 

Plugging (3.59) into (3.53) we then find that 

lab = Cp 2 (g ah + p~ 2 44) (3.61) 

with 

C = e~ 2uh 2-d , p(u) = cosh d ( Uh ~ u 1 . (3.62) 

This is precisely the black brane metric and £ a is the null vector of the horizon hypersurface. 
Consider an arbitrary basis of vectors E m which are orthogonal to £ a , we can expand 

9ab = -l a 4 + oi li E- ia E ]b (3.63) 

Kb = -h 0 £ a h + cr 1} E ia E ib (3.64) 

then equation (3.59) implies that 

cr lk Q!kj = cbi (3.65) 
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where ay is the inverse of a ^ and c is some constant. In other words, regularity condition 
fixes h in terms of g up to two constants c and h 0 . Uh and r can be expressed in terms of c 
and h 0 as 


e du h 


If 2 


e dT = 


4c 1 2 


(3.66) 


c — ho (yfc + -y/ho ) 2 

Given that £ a is the null vector of the horizon and we can choose basis E- w in (3.63) to 
be that in (3.36) (with i index raise and lowered by a) and then a 1 -* of (3.63) then coincides 
that in (3.38), and thus the same notations. Similarly in the horizon limit ho —> 0, and c 1 ' 
of (3.64) is related to cry in (3.36) by raising and lowering using a. 


IV. GENERALIZATION TO HIGHER ORDERS 

In this section we discuss computation of Juv to higher orders in the derivative expansion. 
We first briefly outline the general structure of higher order calculations and then mention 
some results at second order. 

A. Structure of derivative expansions to general orders 

Assuming h , g and X a are slowly varying functions of boundary spacetime variables, we 
can expand 7 ^, the extrinsic curvature /C, and r in the number of boundary derivatives, i.e. 

7 = 7o + 72 H-, /C = /C 0 + /C 2 H-, t = t 0 + t 2 H -. (4.1) 

where 7 0 , /C 0 , r 0 (which we already worked out) contain zero boundary derivatives of 
^,h,J a M = <9 /t X a , whereas 7 2 ,/C 2 ,t 2 contain two boundary derivatives, and so on. One 
can readily see that there is no first order contribution, as the equations for the saddle point 
(2.37) and (2.38) do not have first order terms, and neither does the action (2.32). The final 
hydrodynamical action (1.7) will receive first order contributions as the IR contribution 
will contain first order terms, which will communicate via matching to Tuy at the stretched 
horizons through equations for h\ and h 2 . 
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Let us first look at the dynamical equations (2.37) which under decomposition (3.5) can 
be written as 

(K' + -K 2 ] = -A - - Tr K 2 + d —^ {d) R, (4.2) 

d V 2 / 2 2d y J 

K 7 + KK — TZ— - {d) R (4.3) 

d 

where 77 denotes the matrix of mixed-index Ricci tensor ^R^\. Taking the u derivative 
on (4.2) and using (4.2)-(4.3) we find that 

K" + 3 K'K + (K 2 - d 2 )K = (d) RK + d ~ 2 x ^ R' -— Tr 77/C . (4.4) 

2 (d — 1) d — 1 

Plugging (4.1) into these equations we find at n-th order 

K + :iK oK + (377' + 3 K 2 - d 2 )K n = S n (4.5) 

K + K 0 K n + K n K o = P n (4.6) 

where sources S n and P n contain only quantities which are already solved at lower orders. 
Note that P n is a traceless matrix. Parallel to earlier zeroth order manipulations, the inte¬ 
gration constants in will need to satisfy a constraint from (4.2) 

^ (K' n + K 0 K n ) + Tr K 0 K n = B n (4.7) 

where B n again contains only quantities solved at lower orders. Thus once we have solved 
the nonlinear equations at the zeroth order, higher order corrections can be obtained by 
solving linear equations. In particular, at each order the homogeneous part of the linear 
equations are identical with only difference being the sources. This aspect is very similar 
to the structure of equations in the fluid/gravity approach [10]. For completeness we give 
explicit expressions of various sources S n , P n , B n in Appendix B. 

Similarly at n-th order the r equation of motion (2.28) becomes 

Tr/C 0 /C n - K 0 K n = Y n (4.8) 

where the left hand side should be evaluated at zeroth order solution r 0 and Y n again depends 

only on lower order terms. For example at 2nd order it can be written as 

1 Bt 

Y 2 = --Wr 2 + (V a K 0 - V b !C b 0 a ) ^ . (4.9) 
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Note that r n does not appear in (4.8) as d u (E u u ) 0 = 0. 


B. Non-dissipative action at second order? 


We have carried out the evaluation of Juv to second order. The full results are rather 
complicated and will be presented elsewhere. Here we will only mention results relevant for 
the following question: can we find boundary conditions for h at the horizon which allow 
us to derive a non-dissipative hydrodynamical action to 2nd order in boundary derivatives? 
Mathematically this requires that in taking h to be null, I\jv[h,g, X a ] should have a well- 
define limit and furthermore h will either decouple (as in the zeroth order) or be determined 
in terms of g and X a , with X a unconstrained. There are many reasons not to expect this to 
happen. After all, the holographic system we are working with has a nonzero shear viscosity, 
and things will eventually fall into horizon after waiting long enough time. Nevertheless it is 
instructive to work this out explicitly. Note that ideal fluid action of [23] has been generalized 
to second order in derivatives in [24, 28] based on volume-preserving diffeomorphisms. 

For simplicity we will take g lw = g^ u . We find that in taking the horizon limit A = 
d{uh — r) —* 0, Juv develops various levels of divergences in terms of dependence on A: 


1. The most divergent terms have the form 


A 2 ) 
^uv 


rsj 


trlog 2 a 
A 3 


/ 1 1 
\ log 2 A + log 3 A 



(4.10) 


where cr is traceless part of cdj = a lk a y, and we have suppressed other finite constant 
coefficients. Interestingly all these divergences go away if we impose the regularity 
condition (3.65) at the horizon which is equivalent to 


tr log 2 a = 0 . 


(4.11) 


2. The next divergent terms are of the form 



(4.12) 
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where is the entropy current (3.43). Thus they vanish if we impose the non- 
dissipative condition 

d,r = 0 . (4.13) 

3. Finally, we have the logarithmic divergence of the form 

£ (2) = -^2 3 -%e {2 - d)uh Y, 2 \ogA + O(A 0 ) (4.14) 

Cl 

where 

Z/w = P^ p Puad ip u a) - -^—^d p u p Pp U , Pp U = r) lw + UpU v . (4.15) 

For this divergence to disappear, one then needs 

E 2 = 0 (4.16) 

i.e. the system is shear free. Note that the divergence in (4.14) is proportional to E 2 , 
which is precisely the rate of increase of the horizon area. 4 

If we want to have unconstrained X a , the shear-free condition (4.16) cannot be con¬ 
sistently imposed. Thus it appears not possible to generalize the non-dissipative horizon 
condition to obtain a second order non-dissipative action. As mentioned at the beginning 
of this subsection, this is hardly surprising. In particular, the specific divergence structure 
of (4.14) implies that we must take account of dissipation. 

We should note that in the full Schwinger-Keldysh program (1.7) outlined in the introduc¬ 
tion, there is no need to impose any of the above conditions (4.11), (4.13) and (4.16). The 
divergences will cancel with those from Tip after we do a consistent matching at the stretched 
horizons. Also, the divergences mentioned above are not due to the use of Gaussian normal 
coordinates, which of course become singular themselves at the horizon. Similar divergences 
also occur in Eddington-Finkelstein coordinates. Being off-shell means that there are neces¬ 
sarily both in-falling and out-going modes at the horizon (which will further be magnified 

4 To see this explicitly, one need to study the Raychaudhuri equation associated with the null congruence 
i a on the horizon. In particular, one may need to put on shell the contraction of the Einstein equation 
with i a at this order. See [29] for details. 
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by nonlinear interactions) which will lead to divergences also in the Eddington-Finkelstein 
coordinates. 

Finally, for completeness, let us mention that if we do impose all of (4.11), (4.13) 
and (4.16), we obtain a simple result 

£(2) = 2 l-s e ( 2 ~d) u h 

where 

9 = (5p = P^d u u h , oP = u v d v vP. (4.18) 

The second order Lagrangian (4.17) is subject to the ambiguity in the field redefinition 

(4.19) 

which we fix by using the zeroth order equation of motion 


d — 1 


9 2 -(d- 2)/3 2 - 2 


(4.17) 


P ,ju d,/u h = u v d v vP. 


(4.20) 


to express in terms of aG Eq. (4.17) can then be simplified to 


£( 2 ) 


2 1 ~ : 3e^ 2 ~ d ^ Uh 


~^6 2 - da 2 
d — 1 


(4.21) 


We should emphasize that due to various conditions imposed at the horizon, the off-shell 
nature of the above “action” is not clear at the moment. To have a genuine off-shell action 
for one patch we should first compute the full action for both segments of the Schwinger- 
Keldysh contour, and then integrate out modes of the other patch. At second order this 
“integrating-out” procedure likely does not make sense in the presence of dissipation. Even if 
this procedure makes sense after suppressing dissipation, it is not clear our current procedure 
of imposing regularity and non-dissipative conditions yields that. 


V. CONCLUSION AND DISCUSSIONS 

In this paper we outlined a program to obtain an action principle for dissipative hydrody¬ 
namics from holographic Wilsonian RG, and then developed techniques to compute Juv, as 
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defined in (1.5), at full nonlinear level in the derivative expansion. The “Goldstone” degrees 
of freedom envisioned in [22] arise naturally from gravity path integrals, and the ideal fluid 
action of [23] emerges at zeroth order in derivative expansion when non-dissipative condi¬ 
tion is imposed at the horizon. The volume-preserving diffeomorphisms of [23] appear here 
as a subgroup of horizon diffeomorphisms. We also found that a direct generalization of 
the non-dissipative condition to higher orders does not appear compatible with the action 
principle. 

An immediate generalization of the results here is to compute Tir of (1.3) which will 
enable us to take into account of dissipations. 

In our discussion we have ignored possible corrections from Jacobian in the change of 
variables in going from (2.2) to (2.8), as well as higher order corrections in the saddle 
point approximation of gravity path integrals. Such corrections are suppressed at leading 
order in the large N limit of boundary systems. Nevertheless, they may be important for 
understanding the general structure of the hydro dynamical action, thus it would be good to 
work them out explicitly and explore their physical effects. 

It would be interesting to generalize the results to more general situations, such as charged 
fluids, fluids with more general equations of state (for example [30]), systems with anomalies 
(such as those considered in [31-33]), or higher derivative gravities. 
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Appendix A: Boundary term 


1. Boundary compatible with foliation 

Consider a spacetime M with a boundary DM. Suppose dM is a slice of a foliation of M 
by hypersurfaces We denote the outward normal vector to £„ by n M . The Gauss-Codazzi 
relation gives 

R = + (K 2 - K mn K mn ) - 2 V M (n M V N n N - n N V N n M ) (Al) 

where ^R is the intrinsic scalar curvature of and Kmn is the extrinsic curvature for 
with 

K = V M n M . (A2) 

Now let us consider 

S = [ d d+1 xy/^(R - 2 A) + [ d d x\f^h2K (A3) 

J M J dM 

and apply the Stokes theorem to the last term of (Al) 

— 2 [ d d+1 x Vm(w m V nk n — n N V N n M ) = —2 f d d x V~-hriM(n M 'V n j u n — n JV Vjvn M ) 

J M J dM 

= -2 [ d d x yf—h V M n M (A4) 

J dM 

which then directly cancels the Gibbons-Hawking term. In this case we thus find that 

s= f d d+1 xy/=g [ (d) R+(K 2 - K MN K mn )\ . (A5) 

J M 

2. Boundary incompatible with foliation 

Here we will consider an explicit example with a spacetime metric 

ds 2 = du 2 + 7 ab{u, x a )dx a dx b . (A 6 ) 

We further consider a foliation of the spacetime by hypersurfaces T, u specified by u = const. 
Denote the normal vector field to by n M , which can be written as 

n M = (l, 0 ), n M = ( 1 , 0 ). (A7) 
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(A 8 ) 


The extrinsic curvature for T, u can be written as 

Rah -jd'uHobi A — 7 du^ab • 

The Gauss-Codazzi relation gives 

R = (d) R + (K 2 - K MN K mn ) - 2V M (n M V N n N - n N W N n M ) (A9) 

where ^R is the intrinsic scalar curvature of E u . Now suppose the spacetime region M has 
a boundary dM which does not coincide with one of £ u . More explicitly, we specify dM by 


u = r(x a ) 


(A10) 


for some function r(x a ). The outward normal vector to dM can thus be written as 

1 


£ M = n(l ,-d a r), £ = n(l,-d a r), 

The extrinsic curvature of dM is given by 

1 


n = 


d a T = 7 ab d b T . (All) 


y/1 + d a rd a r ’ 


K I dM — 


d u {ny/^y) - -=d a (ny/^d a T) 


v/=7 V V " 

= n K + n 3 K ab d a rd b r - nD 2 r + n 3 d a Td b rD a D b T 


(A 12 ) 


where we have used (A 8 ) and all indices and covariant derivatives are defined with respect 
to h ab = 7 ab (r(x),x a ). 

The induced metric on dM is given by 


hab h a b “t” d a Tdt>T . 


(A 13) 


Note the relations 

h ab = - h ab + J_ Qa T Qb T ^ ^ = Jfbfo = n 2 ga = n = y/l - d a rd a T 

n 2 

(A14) 

and (A12) can also be written as 

K\q M = nK + -K ab d a rd b T — -D 2 t . (A15) 

n n 

Now combining the boundary term in (Al) and the Gibbons-Hawking term we find that 

S bd = 2 [ d d x\f-hK\ d M-2 [ d d xsfMiK 
J dM J dM 

= 2 f d d x\ /f —Ji-( y (Ji ab (dT) 2 — d a Td b T)K ab — D 2 T ) . (A16) 

Mr n 
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Appendix B: Explicitly expressions of sources 


Here we give explicit expressions of various sources introduced in Sec. IV A 



(Bl) 


i =2 z=2 


1 

/'„ = K n ~ 2 (d>I >n ~ V A ’< K » 


(B2) 


n— 2 n —2 


S„ = Jn -3j2 K 'i K n-i - V !<<!< 


(B3) 



(B4) 


Note that in the last term of (B3), the sum should not include the term with i — j — 0 
(which is denoted using a prime). Also note the relation 



(B5) 
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